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We address the classical capacity of a quantum bosonic memory channel with additive noise, 
subject to an input energy constraint. The memory is modeled by correlated noise emerging from 
a Gauss-Markov process. Under reasonable assumptions, we show that the optimal modulation 
results from a "quantum water-filling" solution above a certain input energy threshold, similar to 
the optimal modulation for parallel classical Gaussian channels. We also derive analytically the 
optimal multimode input state above this threshold, which enables us to compute the capacity of 
this memory channel in the limit of an infinite number of modes. The method can also be applied to 
a more general noise environment which is constructed by a stationary Gauss process. The extension 
of our results to the case of broadband bosonic channels with colored Gaussian noise should also be 
straightforward. 



I. INTRODUCTION 

The growing importance of quantum communication 
motivates a strong research activity on quantum chan- 
nels. The classical capacity of a communication channel 
is defined as the supremum of the rate of classical in- 
formation that can be transmitted via the channel. An 
important question for quantum channels is whether or 
not entangling the channel inputs may improve the trans- 
mission rate. The additivity conjecture proven to be true 
for some memoryless quantum channels (see [l| and refer- 
ences herein) implies that entanglement does not improve 
the transmission rate. Although this conjecture was dis- 
proved recently for another class of channels in Ref . , 
it remains an interesting question in itself to determine 
whether entanglement helps or not to achieve the capac- 
ity. In particular, it was shown for some channels with 
memory, where the environment exhibits a correlation 
be- tween subsequent uses of the channel, that the trans- 
mission rate can be increased by using entangled input 
states. 

In the case of a quantum memory channel with dis- 
crete alphabet, the first studies of the classical capacity 
considered a depolarizing channel [3( and a quasiclassical 
depolarizing channel [H and showed the existence of a 
threshold on the degree of memory above which entan- 
gled signals improve the transmission rate with respect 
to product states. Further studies [!, Q on a qubit chan- 
nel with hnite memory derived bounds on the classical 
and quantum capacities. In Ref. Q, the classical and 
quantum capacities were discussed in a general frame- 
work, where furthermore, a malicious third party may 
have control of the initializing memory. The case of gen- 
eral Pauli channels with memory was studied explicitly 
in Ref. |8j, where it was shown that the optimal states 
are either product states or Bell states separated by a 
memory threshold. The behavior of channels with corre- 
lated error was connected to quantum phase transitions 



in many-body systems in Ref. [9( . The quantum capacity 
of a dephasing channel was computed in Ref. [lij , where 
an enhancement of the capacity with increasing memory 
was proven. For higher dimensions, it has been shown 
that the capacity of qudit channels exhibits the same 
threshold phenomenon as Pauli qubit channels [ill, E3] • 

For a quantum channel with continuous alphabet, it 
was first shown for an additive bosonic channel and a 
lossy bosonic channel, respectively, [H, [3| that in the 
presence of a memory, some degree of entanglement be- 
tween the input states is necessary to achieve the ca- 
pacity, which is in contrast to the behavior reported 
for discrete quantum channels. Indeed, the optimal in- 
put states correspond to Einstein- Podolski- Rosen (EPR) 
states with finite squeezing which in- creases with the 
degree of memory, in contrary to either maximally or 
nonentangled states in the discrete case. A lossy bosonic 
channel with correlated environment and an input energy 
constraint has been studied [HI, [ll| [13] ■ Lower and upper 
bounds for the classical capacity were derived in Ref. |14| 
and the capacity was calculated in Refs. flBL IT5, Il7| ■ 

In this paper, we discuss a bosonic memory channel 
with additive noise modeled by a Gauss-Mark ov p rocess. 
We present an extension of the model of Ref. [l3| , where 
two uses of a bosonic additive channel with correlated 
noise were treated, to the case of n channel uses with 
Gauss-Markov noise. We determine the classical capacity 
of the channel in the asymptotic limit. 

In Sec. [Ill we recall the definition of the classical ca- 
pacity of quantum channels and specify the model of the 
channel under an input energy constraint. Section III is 
devoted to the treatment of a monomodal channel with 
phase-dependent noise where the quantum water-filling 
emerges as optimal solution. In Sec. IV, we introduce 
the method for finding the capacity and apply it to the 
Gauss-Markov memory. Under reasonable assumptions, 
we determine the optimal solution as a global water- 
filling solution and obtain the classical capacity of the 
channel. Section V treats the capacity in several limiting 
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cases: the limit of full correlations, the classical limit, the 
case of a modified thermal noise, and the transition from 
finite to infinite uses. 



II. BOSONIC ADDITIVE CHANNEL 

In order to send classical information through a quan- 
tum channel one chooses an alphabet and associates its 
letters to quantum states p 1 * 1 , where index i labels the 
letters of the alphabet. The input states sent through 
the channel interact with the environment and thus are 
modified at the output. The action of the channel T is a 
completely positive, trace-preserving map acting on the 
input states p\ n : 



T[ P ? 



(1) 



In the messages, each letter appears with a certain prob- 
ability pi so that the overall modulated input state is 
described as a mixture 



P 



E^ 1 



(2) 



By linearity of T, Eq. ([T]) determines as well the action 
of the channel on the overall modulated input 0, that 



is, 



P ; 



(3) 



where we will refer to p as the overall modulated output. 
In order for the overall modulated input to be physical 
it has to obey the energy constraint 



E 



Pi 



Tr^'Va) <n, 



(4) 



where n is the maximum mean photon number per use 
of the channel and will be referred to as "input energy" 
in the following. 

The classical capacity C (T) of the channel T represents 
the supremum on the amount of classical bits which can 
be transmitted per invocation of the channel via quantum 
states in the limit of an infinite number of channel uses. 
This quantity can be calculated with the help of the so- 
called one-shot capacity, defined as in (l8| : 



d(T) 



sup 



X, 



where the Holevo ^-quantity reads 



(5) 



(6) 



with the von Neumann entropy S(p) — — Tr(plogp) 
where log denotes the logarithm to base 2. The supre- 
mum in ([5]) is taken over all ensembles of {pt, pf 1 } of prob- 
ability distributions p t and pure input "letter" states p™ 1 



The term "one-shot" means that only one invocation 
of T is needed to calculate Eq. ([5]). Using this quantity 
the capacity C(T) defined as above may be evaluated in 
the following way. A number n of consecutive uses of the 
channel T can be equivalently considered as one parallel 
n-mode channel T^ n \ which is used only one time. Then 
the capacity C(T) is evaluated in the limit: 



1 



C(T) = lim -Ci{T (n >). 

n—*oc ft 



(7) 



Let us now assume TW to be a n-mode bosonic addi- 
tive channel with memory. In the following, the number 
of modes of this channel corresponds to the number of 
mono-modal channel uses. Each mode j is associated 
with the annihilation and creation operators dj , , re- 
spectively, or equivalently to the quadrature operators 
1i = ( a j + a ])/v / 2, Pj = i{a>j — Oj)/y/2 which obey the 
canonical commutation relation [<7i,£>j] = iSij, where Sij 
denotes the Kronecker-delta. By ordering the quadra- 
tures in a column vector 



R = (qi, ...,q n ;pi, ...,pnf , 



(8) 



we can define the displacement vector m and covariance 
matrix 7 of an n-mode state p as 



where 



m = Tr[pR] 

7 = Tr[(R m)p(R m) T ] - I J, 



1 

-1 



(9) 



(10) 



is the symplectic or commutation matrix with the n X n 
identity matrix 1. In this paper we focus on Gaussian 
states, which are fully characterized by m and 7. Fur- 
thermore, without loss of generality, we set the displace- 
ment of overall modulated states and the means of classi- 
cal Gaussian distributions to zero, because displacements 
do not change the entropy. 

For the bosonic channel, the encoding of classical in- 
formation is made according to a continuous alphabet, 
where the previous discrete letter with index i is replaced 
by the real and imaginary part of a complex number a. 
A message of length n is therefore encoded in a 2n real 
column vector 

a = (SR{a 1 },SR{a 2 },... ) K{a„},S{Q 1 } ) ..,S{Q n }) T . 
Physically, this encoding corresponds to a displacement 
of the n-partite Gaussian input state in the phase space 
by a and is denoted by p 1 ™. The Wigner function of p 1 ™ 
reads 



exp \-\{R- ygq)t 7^ (R - V2a)} 
(2 7 r)™ v /dct ( 7in ) 



(11) 



Throughout this work, we use the recently proven con- 
jecture that a coherent state minimizes the entropy of a 
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mono-modal Gaussian thermal channel [19( . As we show 
in Sees. IIHI and HVl we only need the extension of this 
proof to the case of a mono-modal Gaussian channel with 
anisotropic noise to justify that the optimal n-partite in- 
put state is Gaussian. As a result we only consider Gaus- 
sian distributions of the letters in the messages so that 
the overall modulated input state sent through the chan- 
nel is a Gaussian mixture p m — J d 2n af{a)p™, where 
d 2n a = d3?{ai}d3{ai}...d3fi{a„}d3{a n } with Gaussian 
distribution 



/(«) = 



exp 



-« T 7 m Ld«] 



Tr'^det (7 mo d) ' 



(12) 



centered at zero and characterized by the modulation co- 
variance matrix 7 mo d- 

As we are no longer dealing with probability distribu- 
tions pi but with probability densities f(at), the summa- 
tions in the formulae above are replaced by proper inte- 
grations. The action of T^ n ' on an input state carrying 
a message a reads as in [l3[ 



T {n) {p™}=P° a ut = J d 2 "/3/ cnv (/3) 
x D(p n ) <g> ... <g> £>(&) p™ D*(p\) <8 - ® D\/3 n ), 



(13) 



with d 2n (3 = «{/3i}dS{/3i}...d5R{/3„}d3{/3„}, j3 = 
(fc{/3 1 },...,fc{f3 n },%{l3 1 },..., < S{(3 n }) T and the displace- 
ment operator D((3j) = e^°i - ^* °j '. The displacement 
is applied according to the Gaussian distribution of the 
environment 



/env(/3) 



exp [-/3 T 7^ (3] 

TT n y/det (7env) 



(14) 



with (classical) covariance matrix 7 env . If this matrix is 
not diagonal, then the environment introduces correla- 
tions between the successive uses of the channel. These 
correlations model the memory of the channel. 

Since we centered the distributions of the environment 
and modulation as well as the Wigner function of the 
input state around zero, the covariance matrices of the 
state carrying the message a at the output of the channel 
p° ut and of the overall modulated output state p, read, 
respectively, 



7out Tin T" 7env 
7 — 7out T" 7mod - 



The one-shot capacity of such a system is 
d(T("))= sup X n, 

7in,7mod 



(15) 



(16) 



with the Holevo %-quantity which then reduces to 

Xn=S(p)-S(p°» t ). (17) 



In the case of a Gaussian state p, the von Neumann en- 
tropy can be expressed in terms of the symplectic eigen- 
values i>j of its covariance matrix: 



S(P) = J> 

g(x) = 



(18) 



(x + 1) log(x + 1) — xlogcc ,x>0 
, x = 0. 



The energy (or mean photon number) constraint (UJ) is 
given by 

^ (Tr ( 7in ) + Tr ( 7mod )) -\=n. (19) 

We remark that in Eqs. (|17H19[) the quantities depend 
only on the covariance matrices and not on a particular 
message a. Therefore, we can fully discuss the action of 
the channel solely in terms of covariance matrices. 



III. MONO-MODAL PHASE DEPENDENT 
CHANNEL 

First we study the mono-modal channel T, where the 
input state undergoes a phase-dependent noise with co- 
variance matrix 



7env 



V 9) 

fenv 

, o 



o 

/env , 



~(9) 
/ env 



^7, 



(p) 

env ' 



(20) 



The optimal solution of a channel with such a noise has 
not yet been fully studied. In Ref. [2(| the capacity 
of a multi access channel with such a noise was studied 
but the optimal input and modulation were not discussed 
in detail. It was shown [l9} for a mono-modal thermal 
channel that a coherent input state which is modulated 
according to a Gaussian modulation achieves the capac- 
ity. We conjecture that this proof can be extended to the 
case of the noise given by Eq. (|2H)) . i.e. that a Gaussian 
input state and a Gaussian modulation remain optimal. 

The general input and modulation covariance matrices 
read 



,-v. I 'in /in 

T in - (w) w 

\ /in 'in , 



/ (?) <«0\ 

, — /mod 'mod ('OlA 
\ 'mod 'mod / 



We now determine the optimal 7i n , 7 mo d by solving 
equation Eq. (Til))) with energy constraint p^|) and the 
requirement to have a pure input, i.e. 



det 7i, 



(22) 



With these constraints we can write out the total La- 
grangian of the system, i.e. 



4 



C = 



C(7^ ) +7^ d +7^ ) +7^ d )-r(7^ ) ), 



(23) 



r 



'mod 



Jq) 




Quadratures 



FIG. 1: (Color online) Quantum water-filling solution for the 
mono-modal channel: Optimal eigenvalues for both quadra- 
tures. 



where £ and r are Lagrangian multipliers. We arrive at 
a system of eight equations where the solution reads as 
follows. First, we find that 



7 



iip) 



(IP) 

' mod 



o. 



(24) 



i.e. all three covariance matrices are diagonal in the same 
basis. Secondly, we obtain 



Ji) + 7 (?) + 7 («) 

/in 1 'cnv 1 /mod 



'in 1 'env 



(p) 

Tmod ' 



(25) 

which can be regarded as a "quantum water-filling" [28j 
solution and implies that the optimal output state 7 is 
a thermal state (see Fig. [T]). This repartition of input 
energy for the preparation of the quantum state and the 
classical modulation generalizes the situation of parallel 
classical channels with a joint energy constraint. Here, 
we have two channels, one for the q quadrature and one 
for the p quadrature. The eigenvalues of the covariance 
matrix of the modulated output state read 



7 



(?) 



7 



(p) 



/cnv 



/cnv 



(26) 



which simply reflects the total energy of the system. 
Thirdly, we obtain the optimal degree of squeezing of 
the input state: 



7 



(?) 

in 

» 



/cnv 
/cnv 



(27) 



Thus, the squeezing at the input has to match exactly 
the phase dependence of the noise. 

Since all matrices are diagonal we conclude that 
the variables 7/* 7^,7^ od >7raod are tne eigenvalues of 
7in,7mod- Using the condition that the input has to be 
pure (|22|) we conclude from Eq. ([27]) that 



(?) = I 

7 in 2 ^ 



1 



7env (p) _ _ 

rM ' 2 

/env 



\ 



/cnv 
/env 



Then, the modulation eigenvalues simply read 



' mod 



Tout 



' mod 



(p) 

Ai-7out> 



(28) 



(29) 



where 7 ou f } 



7i 



(?>p) 



(?>p) 

/cnv 



and where the quantum 
water-filling level [i = 7^ = 7^ is determined by 
Eq. ([2H]). Equations (gHJ and ([Ml) are depicted in Fig. [f[ 
Note that both modulation eigenvalues (|2"9"]l have to be 
non-negative to be physical. This condition together with 
Eq. ([28]) result in the threshold on the input energy 



n > nthr = 



1 



max { 7cnv,7cnv} 

.(«) >)| 
env 5 /cnv J 



^ I \ min {7, 



+ 17, 



(9) 

cnv 



(30) 

below which the obtained solutions do not hold. 

Knowing all optimal eigenvalues of the system, we 
can easily determine the one-shot capacity defined by 

Eq. com, 



/cnv 



/env 



^-y 1 , /env 1 /env 1 / / 101 \p\ \ ^ 

^i = 9 I n H I -5 I V 7env 7env ) • n > 7l t hr- 

(31) 

One notices that in the first term the algebraic mean 
of the noise (7env + 7env)/2 and in the second term the 
geometric mean [7cnv7onl] -1 ^ 2 appears. Comparing C\ 
with the definition of the Holevo quantity ^ one can 
say that the first term corresponds to the entropy of a 
thermal state with a mean photon number identical to 
the mean energy of the system, and, roughly speaking, 
the second term corresponds to the average of the en- 
tropies of two thermal states with means 7e~nv,7cnv. In 
the case 7!^ 



/cnv 



N we recover the capacity of a 
mono-modal thermal channel C = g(n + N) — g(N). 
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IV. GAUSS-MARKOV CHANNEL 

Now we proceed with the central part of this paper: 
the treatment of the multimode channel as introduced in 
Eq. ([13")) with a Gauss-Markov noise. 



A. The noise model 

Assume a Gaussian distributed real random vector Z 
that is generated by a Markov process: 



>Zi-i+W u <</><!, 



(32) 



where Wi are Gaussian distributed and identically and 
independently distributed (i.d.d.). We set without loss 
of generality the expectation values E(Zi) = E(Wj) = 0. 
We request in addition that the variance Yax(Zi) = N, 
which leads to Var (Wi) — (1 — <fi 2 )N. Then we conclude 
that the matrix elements of the covariance matrix M for 
the Gauss-Markov process (|3"2"]) read 

MyO) ee Cov(Z 4 , Zj) = N0 lt ~ jl , <cf>< 1, (33) 

where N £ R denotes the variance and (j) is the near- 
est neighbor correlation. Equation (f3"3"]l is a symmetric 
Toeplitz matrix which is completely defined by its diag- 
onal elements £^ M ^ = Nc/)\ k \ — r^O (see Appendix |A]) . 

In the following, we treat the channel given by Eq. JT - 
with a classical noise covariance matrix 



7env 



_ (M{4>) 



M(-(j>) 



(34) 



where M((f>) reads as in Eq. ((33|) . In other words, the 
noise of the q (p) quadratures results from a Gauss- 
Markov process with (anti-)correlated nearest neighbors. 



B. Solution 

In order to determine the capacity of the memory chan- 
nel (|13p . one needs to find the optimal input and mod- 
ulation covariance matrices which obey the energy con- 
straint (TH)]) . One can treat equally this problem in the 
basis where the channel becomes memoryless provided 
that one can rotate the noise covariance matrix by a pas- 
sive symplectic transformation into the basis where it is 
diagonal. This is possible, because a passive symplectic 
transformation does neither change the entropy, nor the 
energy constraint of the system. For the class of addi- 
tive bosonic channels, where the noise covariance matrix 
7env does not contain any cross correlations between the 
quadratures, i.e. has the form 



7e 



y«) o 

-v (p) 

. u fenv 



(35) 



one can show, that the transformation that diagonalizes 
the noise is symplectic only if the two block matrices 7env 
and 7envj respectively, commute (see appendix |A"| . 

Applied to the Gauss-Markov noise 7 env as defined in 
Eq. ((Ml), this condition requires that [M (</>), M(-<j>)] = 
0. However, for finite n one can show that this does not 
hold and, as a consequence, 7 e nv cannot be diagonalized 
by a passive symplectic transformation. In order to over- 
come this problem, we introduce another noise covariance 
matrix ~f' em that does have commuting block matrices 
(see (|35p ) and, furthermore, an eigenvalue spectrum that 
asymptotically converges to the eigenvalue spectrum of 
7env The capacity of this modified channel system and 
of the original one will be identical, since this quantity is 
found in the limit of an infinite number of channel uses. 

We introduce the noise covariance matrix 



7e 



~ I m( c )h 



where M^ c > ((f)) is a circulant symmetric matrix 



< \i —j\ < k 
n < \i — j\ < n 



1. 



(36) 



(37) 



with k = (n— 1)/2 for odd n and k = n/2 for even n, and 
i,j = 1, n. It is known that all circulant symmetric 
matrices commute, and furthermore, that the spectrum 
of a symmetric Toeplitz matrix and a circulant symmetric 
matrix generated by the same diagonals asymptotically 
coincides (see appendix [Al for further details). 

At this point, we base the following calculation on the 
conjecture that the optimal input covariance matrix "f in 
diagonalizes in the same basis as the noise covariance 
matrix 'y' env . We rotate the matrix 74 nv into the basis 
where it is diagonal and therefore treat a set of n un- 
corrected channels. Then, for each channel the optimal 
input takes the form derived in Eq. (|28p. provided the to- 
tal input energy suffices to exceed the threshold for each 
channel. We note in addition, that by Eq. (jT7j) and the 
subadditivity of the entropy 



S(p) < Sfa) + S(p 2 ) 



where 



Pk 



Tr- 



Pll — lPk-llPk+H 



(38) 



(39) 



the optimal modulation covariance matrix ^' moA is also 
diagonal in this basis. 

We now take the limit n — > oo. In appendix [Bj we 
determine the spectrum of M(<fi) in this limit, which is 
identical to the asymptotic spectrum of M^°'((j)). Thus, 
as the spectra of the two noise matrices 7 en v and ~/' ellv are 
now identical, the optimal input and modulation covari- 
ance matrices y ln and j' in , (7^ od and j' mod ) coincide. 

In this limit the spectrum of 7 cnv as introduced in 
Eq. ((34]) becomes continuous and reads for the two 
quadrature blocks 



7 (<r) 

' env 



(x) 



l env 



(7T- 



A( M '( 3 



(40) 
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where \( M ^(x) is the spectrum of M{4>) given by 
Eq. (TBTj), i.e. 



\( M \x) = N 



l + 2 -2</>cos(ir)' 



(41) 



where x G [0, ir] is the spectral parameter. 

For each individual channel, we know from Eq. 
that if its allocated input energy is sufficient then the 
optimal overall output state is a thermal state, i.e. 



As in the mono-modal case each input has to match the 
anisotropy of the corresponding noise, i.e. 



(</) 



(p), 



» 



modulation energy 
per mode (1-r|)n 



Spectral parameter x 



(,•) 



(43) 



In order to check whether the optimal input state is 
entangled in the original basis (i.e. the basis where the 
noise covariance matrix is introduced in (|34p h we use 
the transformation Q given by (|A4[) for each quadrature 
block, to rotate the covariance matrix of the state back. 
In this basis, we conclude that the first mode is in a 
thermal state for correlations <fi > (a coherent state for 
4> = 0), with variances 



FIG. 2: (Color online) Global quantum water-filling solu- 
tion: Output spectra ~fiu t (x) and ^^(x) shown by the solid 
decreasing curve and dashed increasing curve vs. spectral 
parameter x. The solid bar represents the global quantum 
water-filling level jU g i; the light gray area shows the energy 
used for modulation in q. 



where J^t\x) = J^' p \x) + Tenv* (x) is the (quantum) 
output spectrum and figi is now a global quantum water- 
filling level over all channels, such that 



'in,l 




2(f) cos (x) 



(44) 



Since the purity of the overall input state is untouched by 
the rotation, we conclude that the first mode is entangled 
with the rest of the modes. A full proof of entanglement 
could possibly be obtained by the same method applied 
to all modes. 

We define the fraction of input energy used to prepare 
the squeezed (or entangled) input by 



dx^{x) = {l-ri)n. 



(47) 



The solution is graphically represented in Fig. [5J One 
concludes that the input energy threshold which deter- 
mines whether the modulation (|46p is positive over the 
whole domain x € [0, n] is given by 



(1-Tj)n> 7 &>(0) 



dx^ t (x) 



(48) 



(45) 



\0 



where n is the mean maximum photon number constraint 
defined in Eq. . 

From classical information theory, we know that the 
amount of classical information sent through a system 
of parallel channels is maximized by a water-filling solu- 
tion [HI . By Eq. (|42|) we confirm that this holds for the 
optimal modulation spectrum, i.e. 



Using the definition in Eq. (l4"S"l) and the property of the 
environment 



- f dxj^(x)=N, 



we define 



«thr 



(49) 



(50) 



i.e. the minimum input energy needed to modulate the 
whole channel spectrum. 

If we assume that for given <j>, N 



7mod ( x ) = Mgi 



(46) 



n > n t hr(0, N), 



(51) 
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then (see FigHJ, using Eq. (J45J) and Eq. (gU), we deter- 
mine the global quantum water-filling as 



Mgi = (1 - rj)n + i / dx ^ t {x) 
o 



(52) 



which corresponds to the overall modulated output vari- 
ance of each of the channels, i.e. 



1 



^)(x)=^»(x) = » sl = n + N+~, (53) 

where j( q ' p \x) correspond to the spectra of the covari- 
ance matrix 7 as defined in Eq. (fT5|) . With the definition 
of the symplectic spectra for the quantum output, the 
total output and the noise in the asymptotic limit, that 
is 



V(x) = ^^\x)^p\x) (54) 

V cm {x) = \J 7cnl(x) 7enl(^), 

we are now ready to determine the capacity of the channel 
above threshold for an infinite number of uses: we take 
the limit of Eq. (fTB"]) where the sum tends to an integral 
over the whole domain, i.e. 



C = lim 



1 



n— >oo n -y. 



SUp Xn 



7T 

~ j dx lg (v(x) - - g (v out (x) - j 


g {n + N)-- [dx g (vS^T&kT) • (55) 



g(n + N) / dx g (u env (x)), 



(56) 



A. Classical limit 

In the classical limit we increase the input energy and 
the noise while keeping the signal to noise ratio constant. 
From the limiting behavior of g(x): 



lim [g(x) — log(ir)] = 0, 



(57) 



we conclude that we have to replace the ^(x)-functions 
in Eq. (|56|) by logarithms. The integral term can be 
simplified with the help of Ref. [22| 



7T 

^Jdx log ^ 7 &l(a:)7&l(30) = log(JV(l - 2 )). 



Thus, we conclude that 



(58) 



lim 

\N—*oo, n/N- 



C = C ci = log 



I -ft 2 



n 



Since g(x) < log (x) we reach this limit from below and 
hence, for finite n, N, the classical capacity of the bosonic 
channel is always smaller than the capacity of the clas- 
sical Gaussian additive channel. This is the expected 
result because for the bosonic channel (for <j> ^ 0) a cer- 
tain amount of energy is needed to prepare the input 
squeezed state, which is not the case for the classical 
channel. In the classical limit, the relative fraction of 
energy used for squeezing vanishes to zero and there- 
fore the capacities coincide. In Fig. [3] we plotted the 
capacity C and relative fraction g versus the noise pa- 
rameter N. For each <j> we fixed the signal to noise ratio 
n/N = nxhr{(t>,N')/N' ', with N' = 1 which guarantees 
to be above threshold ViV > 1. The capacities all tend 
to their classical limit (I59|) . while all fractions g decrease 
towards zero. Interestingly, for fixed signal to noise ra- 
tio, we find that g{4> = 0.7) > g((j) = 0.9). This is due 
to the fact that n t hr(0 = 0.7, 1) < n t h r (4> = 0.9, 1). Al- 
though we need more input energy to match the higher 
correlation its cost with respect to the total energy n is 
lower. 



where v cm (x) is the asymptotic symplectic spectrum of 
7env We conclude that the capacity is given by the dif- 
ference of the entropy of a thermal state with mean pho- 
tons n + N and, roughly speaking, the mean entropy of 
the environment. This is a generalization of the capacity 
found for the mono-modal channel found in Eq. (|3ip . 



V. LIMITING CASES 

In this section we analyze the capacity (|56p in several 
limiting cases, such as the case of the classical limit, the 
limit of full correlations, the case of a symmetric noise 
model and the transition from a finite to an infinite num- 
ber of channel uses. 



B. Full correlations 

We now investigate the limit of full correlations, i.e. 
(j) — > 1. Since the first term in Eq. (|56[) is independent of 
</> we take only the limit of the integral term (where we 
drop the constant N inside the eigenvalue functions for 
now). For the integrand we find that 



lim q 



1 



0, < X < IT. 



(60) 

At the borders x = {0, it} the integrand is equal to one 
for arbitrary (j). As this contribution is however infinitesi- 
mally small it is easy to show that the integral in Eq. 
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FIG. 3: (Color online) Fraction of input energy used for 
squeezing rj and Capacity C in bits per use vs. noise parame- 
ter N. The increasing curves correspond to the capacities, the 
decreasing to r\. The gray horizontal lines correspond to the 
classical limit C c \ of each capacity. For all graphs the dotted, 
dashed and solid line correspond to <j> = 0.4, 0.7, 0.9. For each 
<j> we set ViV > 1: n/N = n thr (0, N')/N' with N' = 1 . 



vanishes and we conclude that 



lim C = g(n + N). 

A— >1 



(61) 



This result is in contrast to the result in the classical 
limit (l59l). for which 



lim C c \ 



(62) 



Note that Eq. (|6ip is valid only when the waterfilling con- 
dition (|5D|) is satisfied. However, as <f> gets close to one 
while n is fixed, at some point this condition will be vio- 
lated. In order to satisfy this condition with increasing <fi 
one has to increase n according to (|50|) , so that ([BTjl be- 
comes correct only as an asymptotic limit, with the right 
hand side diverging as goes to one. Nevertheless, we 
conjecture that for fixed n the capacity of the quantum 
channel will be finite in the limit (j> — > 1, because even in 
the absence of noise the capacity is finite as it is equal 
to g(n). Hence, unlike in the classical case, for our quan- 
tum channel there is no diverging benefit from increasing 
correlations. This conjecture will be studied further in a 
forthcoming paper. 



C. Symmetric correlations 

In this subsection, we consider a modified, symmetric 
environment with same correlations in both quadratures: 



7e 



_ (M{<t>) 
- 1 M(4>) 



(63) 



As both quadrature blocks have now identical spectra, we 
treat a set of independent thermal channels, when 7 en v,s 
is diagonalized. We recover from Eq. (|43ll immediately 
that in the asymptotic limit the optimal input spectra 
read 



\ 



-Y (9 ' P) (r) 
fcnv.s V x } 

(p,q) 

/cnv.s 



(,:) 



1 

2' 



(64) 



i.e. the overall optimal input state is a set of coherent 
states and entanglement does not improve the transmis- 
sion rate. This is in agreement with previous investi- 
gations of the two- mode model discussed in [13]. Re- 
cently, the same observation was made independently in 
Ref. Hi. 



D. From finite to infinite number of uses 



If one applies the method introduced in section IIVBI 
to a noise covariance matrix, where the two quadrature 
block matrices (see (J2S])) commute for finite n, then the 
optimal transmission rate reads as in (|56p , where the in- 
tegral is replaced by a sum term normalized by n: 



R (n) = 9(r 



n ^— ' 

fc=i 



l {q \l {p \ 

'cnv.fc 'cnv.A; 



(65) 



where 7^ k , 7^ k denote the eigenvalue spectra for the 
given environment. 

Although, for finite n the two block matrices M((j>) 
and M{—<f>) of the introduced Gauss-Markov noise (|3~4"]) 
do not commute, we investigate as it converges with 
increasing channel uses n to the capacity (|56|) . where 

7cnt,fc'7™t,fc are the spectra of M ((/)), M (-</>) and nu- 
merically obtained. In Fig. [4] we plot for fixed N, 
various correlation strength <f> and a fixed input energy n, 
which is above threshold for the strongest (j>. In addition 
we denoted by gray bars for each 4> the asymptotic capac- 
ity given by Eq. (|56|) . In the plotted region, we observe 
that R^ indeed converges to the capacity C. 



VI. CONCLUSIONS 

The classical capacity of a multi-mode channel with 
Gauss-Markov noise was found under certain assump- 
tions above a certain input energy threshold. By di- 
agonalizing the noise covariance matrix one was led to 
treat first the capacity of a one mode phase-dependent 
channel. 

In the one-mode case we have shown that above the 
threshold the optimal quantum input state is a squeezed 
state where the squeezing matches the anisotropy of the 
noise. Furthermore, we found that the optimal overall 
modulated output state is a thermal state and therefore 
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FIG. 4: Function and capacity C in bits per use vs. 

number of channel uses n. The dotted, dashed, dashed-dotted 
and solid curve correspond to R (n) with = 0, 0.4, 0.55, 0.7. 
For each <f> the corresponding capacity is shown by a gray 
horizontal line. For all plots we took N = 1 and n = 7.5 > 
n thl (<j> = 0.7, N = 1). 



appendix [X] for more details). This means that for any 
pair of Toeplitz matrices T, T" replacing M ((/)), M(—<p) 
in the capacity (above threshold) reads as in (|5rj|) 
where the noise spectra are replaced by the spectra of 
T,T'. 
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APPENDIX A: TOEPLITZ AND CIRCULANT 
SYMMETRIC MATRICES 



the classical modulation is determined by a "quantum 
water-filling" level. 

For the multi-mode channel with Gauss-Markov noise 
the optimal input and modulation were discussed in the 
asymptotic limit. Above an energy threshold, the op- 
timal input eigenvalue of each channel is similar to the 
mono-modal solution determined by the anisotropy of the 
noise of the channel. The optimal modulation was found 
by a global quantum water-filling, which led us to the 
conclusion that the overall modulated output state is the 
same thermal state for all channels or equivalently all 
channel uses. When we rotated the covariance matrix 
of the overall input state back to the original basis, we 
confirmed that the first mode is in a thermal state and 
hence the total state is entangled. 

Finally, several limit cases were discussed. First, the 
classical limit of the channel was discussed and we showed 
that the quantum expression tends to the classical ex- 
pression. Secondly, we argued that in the limit of full 
correlations the capacity of the quantum model stays fi- 
nite whereas it is diverging for its classical counterpart. 
In the case of a symmetric noise model, we recovered that 
the optimal input is given by a set of coherent states. At 
last we depicted the asymptotic behavior of a function 
which represents the optimal transmission rate for finite 
uses for commuting quadrature block matrices, and ob- 
serve that it tends to the capacity for infinite channel 
uses. 

In addition we remark that the proposed solution 
method can be applied more generally to an environment 
which is constructed by a stationary Gauss processes, 
because these processes are characterized by symmetric 
Toeplitz matrices which commute asymptotically (see the 



In this section we specify the symplectic transforma- 
tion that diagonalizes the noise 7g nv (|36|) for finite di- 
mension n and 7 en v (|34[) for n tending to infinity. 

A quadratic matrix T with dimension n x n is called 
Toeplitz Q if 



T 



t 



(T) 



(Al) 



(T) 

where t\_'- is a series of real numbers and i, j = 1, . 



1-3 



The matrix in Eq. (|A1[) belongs to the Wiener class, if 
i r r\ 

{t k } is absolutely convergent, that is, 
J2T=-oo l^^l < 00 • Then the Fourier series 



{T) ikx 



x e [0, 2tt1 



(A2) 



k=— oo 



exists and is Riemann integrable. 

A Toeplitz matrix of dimension n is called circu- 
lant symmetric [25j if it has the form 



T (C) = 



(to ti ti ■ ■ ■ ti t{\ 
t\ to t\ ■ ■ ■ £3 t-x 
ti t\ to ■ ■ ■ ti t$ 



(A3) 



Vl t2 t^ ■ ■ ■ ti to J 

The unitary transformation Q, that diagonalizes ' 
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with odd dimension n, reads 



V n 



V 2 


1 

To 


i 




1 


cos(^) 


cos(^) • 


, 2-Wn-l) \ 
COS ( i '- ) 





sin(^) 


sin (^) ■ 


sin ( — 1 ) 


1 


cos(^) 


cos(^) • 


COS ( i 1 ) 


V o 




sin(^if) • 





(A4) 

For even n the same pattern holds, except for a row 
n -1 ' 2 (l, — 1, 1, — 1) at i — n/2. The diagonal matrix 
Q T T^ c ' Q converges with increasing n to D, where 



D = diag(rfi,d 2 , ...,d v 



where for odd n 



l 2j 



di 
hj+i 



MO) 



(A5) 



(A6) 



d n = /t(tt) for even n, where j — 1, 2, (n — l)/2 and 
/t(x) as defined in (IA2|) . In addition, if a Tocplitz ma- 
trix T belongs to the Wiener class, then Q T T Q also 
converges to D. This means that T^ c > and T can be 
asymptotically diagonalized in the same basis, and hence 
that all Toeplitz matrices (which belong to the Wiener 
class) asymptotically commute. 

Now we show, that the transformation, that diagonal- 
izes 7en V as defined in (|36|) is indeed symplectic. A trans- 
formation S is called symplectic if 

SJS T = J, (A7) 

For 

(A8) 
(A9) 



where J is the commutation matrix defined in (|T0 
a covariance matrix 7 of the shape 



7 



l q 

7 P r 



a rotation R that diagonalizes 7 is of the shape 



1 Ur 



Therefore, in order for R to be a symplectic transforma- 
tion we find the requirement that 



ui u p = u q u; 



(A10) 



where U q ,U p are orthogonal transformations. For the 
noise (|36[) . both quadrature blocks 7q, 7 P are diagonalized 
by the same orthogonal transformation U q = U p = Q 
and thus R corresponds to a passive symplectic transfor- 
mation. 



APPENDIX B: SPECTRUM OF A 
GAUSS-MARKOV COVARIANCE MATRIX 

Our goal in the following is to determine the spectrum 
of M(4>) defined in Eq. (|33[) . As this matrix generates 
a Markov process, one intuitively expects Af _1 (0) to be 
three diagonal, such that only nearest neighbor terms 
appear in the exponential of the resulting Gaussian dis- 
tribution. Therefore we consider a matrix 



/ 1 



V 



1 1 



N 1 -c/) 2 







1 + 02 
1 









1+0 2 
1 



1+v 



where ^ N £ R, < \4>\ < 1. We observe that 



K — --OQ 

and therefore obtain the Fourier series 



< 00, 



7 

(Bl) 
(B2) 



fv(x) = 



Nl 



1 



1 



cos(x) , (B3) 



where the right hand side of the latter was obtained by 
using (|A2[) . From [24[ we state that if the spectrum of 
a Toeplitz matrix T is strictly positive, then the inverse 
T~ 1 is asymptotically Toeplitz, with diagonals 



"T'-s'"' 



-ik'x 



Inserting f|B3|) in the latter leads to 

or equivalently, we found that 

lim V M~ 1 (4>). 



(B4) 



(B5) 



(B6) 



Hence, in the limit of infinite number of rows and columns 
we found the inverse of Mj6). We determine the spec- 
trum of V by using Ref . [26| , take its inverse and there- 
fore receive the spectrum of M{<p) in the limit of infinite 
channel uses: 



1 



1 + 4> 2 - 2^cos(a;) : 
with the spectral parameter x € [0,7r]. 



(B7) 
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